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1 . INTRODUCTICM 


This document describes the implementation details of that part 
of the CAMS Development Dot Data System whicl. has been put under 
formal configuration control. 

This consists of a data base, eight programs which can be used 
to maintain the data base, and fifteen subroutines which can be 
used in Fortran progreuns to process the data in various ways. 
This system supports the task of developing products to aid the 
analyst interpreter in identifying grid dots. Since this is a 
development system, the ability to quickly and easily change the 
application software is important. The programs described here 
contribute to this end. 
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2 . APPLICABLE DOCUMENTS 


TIRE 79-0004, "File Storage Problem", March 1979. 
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3. SYSTEM DESCRIPTION 


3.1 GENERAL 

The CAMS Development Dot Data System (CDDDS) has 3 parts. The 
first is the data file and its index. The data file consists or 
records containing 256 Integer *2 words, usually spectral values 
for one channel for 209 dots or labels for 209 dots. The index file 
contains segment pointers for the data file. The second part of 
the system is a set of programs which can bi used to maintain 
the data files. Third is a collection of subroutines which can 
be called by programs to use the stored data for exploratorv and 
development tasks. 

3.2 FILE DESCRIPTION 

The data files for the CDDDS consist of 512 byte records, a size 
chosen to minimize I/O processing time (Fortran IV - Plus User's 
Guide pg 4-9) . In order to avoid confusion with other files, 
these files were given the type .ACA, by A1 C. Anderson, their 
designer. There are two different types of .ACA files, index 
files and data files. 

3.2.1 THE INDEX FILE 

An index file consists of duplicate record pairs. Each record is 
organized as 256 Integer *2 words. The first record pair, record 
1 and 2 contains only 2 words: 


Word 

1 - The 

number 

of 

records 

in 

the 

index 

file 

Word 

2 - The 

number 

of 

records 

in 

the 

data 

file 


All other words are zero. 

The subsequent identical record pairs contain segment identifiers 
in odd words and pointers in even words: 
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Word 

2n-l 

The segment number 
stored in the data 

for a segment which 
file (or zero) . 

has been 

Word 

2n 

A record number in 

the data file which 

contains data 


for the segment number in word 2n-l (or zero) . 

One of these files may have pointers for 128 segments. An 
index file typically has 10 or fewer records. 

3.:^. 2 THE DATA FILE 

Each record of the data file contains 256 Integer *2 words. 
209 of the words typically contain spectral values fee one 
channel for 209 dots or labels for the 209 d..ts. The other 
47 words contain identifying information, data base pointers 
and ancillary information. 

The input sources do not always contain the ancillary informa 
tion, so the ancillary information may be in error. 

Data file record format. 

Word 

1 Segment identification 

2 Acquisition year 

3 Acquisition day 

4 0 

5-7 Record type code 

8 0 

9-217 Dot specific data 
218 0 

219 Sun angle 

220 Soil line greenness estimate 

221 Satellite identification 

222 Latitude degrees 

223 Latitude min. 

224 Longitude degrees 

225 Longitude min. 
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6-253 

0 

254 

Delete flag 

255 

Pointer -back 

256 

Pointer-forward 


1 


3.3 PROGRTIMS FOR PILE MAINTENANCE 
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3.3.1 LOAD LABELS PROM CARD (CDLABLD) 

3. 3.1.1 Linkage 

This program calls GETNAM entry GETCRD and PUT. 

3. 3.1.2 Interface 

The interface with GETNAM is through the calling arguments the 
interface with PUT is through the calling arguments and the 
common blocks ANCIL and FILNAM. 

3. 3. 1. 3 Input 

Data is read from cards in the default file F0ROO1.DAT. 

For each segment to be loaded there must be a HEADER card with 
format ; 

COL 1-6 HEADER 

COL 10-13 Segment number 

COL 15-20 Code name, for example TRTH78 

HEADER 1663 TRTH78 

The dot labels may be entered in two different formats on DOT 
or LINE cards. 


DOT card format (like ERIPS format) ; 

COL 1-4 DOTb (b is blank) 

COL 6 Label 1 

COL 8 Label 2 

COL 10-80 Dot numbers separated by one or more blanks. 


For example, the card DOT W 1 37 123 201 assigns label W1 

to dots 37, 123, and 201. 
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LINE card format: 


COL 

1-4 

LINE 



COL 

10-11 

Line number ( 1 , 2 . . . , 11 ) 

COL 

13-14 

Label 

for 

first dot of line 

COL 

16-17 

Label 

for 

second dot of line 

COL 

67-68 

Label 

for 

19th dot of line 


For example, the card 
LINE ^ W1 A3. 

Assigns label Wl to dot number 96 or line 6 col 1. and label A3 
to dot 97. 

The last data card should have STOP in col 1-4. If this is 
missing, the statement CDUkBLD STOP BAD DATA? will print on 
the operators terminal at the end of the job. 

3. 3. 1. 4 Output 

This program creates a report (fig. 3. 3. 1.4) and writes new 
records into the ACA data tie. 

3. 3. 1.5 Storage 
955 Bytes 

3. 3. 1.6 Description 

This program reads cards for one segment, loading the labels 
into an array. When it encounters the next HEADER card or a 
STOP card it calls PUT to load the data into the ACA data base. 
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3.3.2 LOAD SPECTRAL VALUES PROM CARD (CDLOAD) 

3. 3.2.1 Linkage 

This program calls GETNAM entry GETCRD, PUT and the private 
subroutine HEX . 

3 . 3 . 2 . 2 Interface 

The interface with GETNAM is through the calling arguments. 

The interface with PUT is through the calling arguments and 
through the common blocks ANCIL and FILNAM. 

The interface with HEX is through the calling arguments. 

3. 3. 2 . 3 Input 

The input cards are read from the default ^ile F0Rf/Ol.DAT. The 
first card has the code identifier (e.g SWCHl*) in col 1-6. 
Following is one or more decks generated by the Spectral 
Unload Utility on ERIPS. 

3. 3. 2. 4 Output 

The output is a line printer report (fig 3. 3. 2. 4) and new records 
in the ACA data file. 

3.3. 2.5 Storage 
2599 Bytes 

3. 3. 2.6 Description 

The code name is read. Then spectral values are read, translated 
by HEX from hexidecimal alphanumeric to decimal integer and 
loaded into 4 arrays. When the arrays are full, they are loaded 
into the ACA data file by PUT and a new set of vlaues is read. 
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3.3.3 DUPLICATE FLAGGING (PURGE) 

3. 3. 3.1 Linkages 

Purge calls GETNT^ entry GETCRD. 

3 . 3 . 3 . 2 Interface 

The interface with GETNAM is through the calling arguments. 

3. 3. 3. 3 Input 

The ACA index and data files as named in FN.DAT are read. 

3. 3. 3. 4 Output 

Duplicate records in the data file are marked and a Line Printer 
report (fig. 3. 3. 3. 4) is generated. 

3. 3. 3. 5 Storage 
1795 Bytes 

3. 3. 3. 6 Description 

The records for each segment are searched for duplicates. \AThen 
a duplicate is found, byte number 254 is set to .TRUE. Otherv;ise 
byte 254 is set to .FALSE. Other programs in the system ignore 
records marked .TRUE. 
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3.3.4 INTERACTIVE QUERY (QUERY) 

3. 3.4 .1 Linkage 

This program calls GETNAM. 

3. 3.4.2 Interface 

The interface with GETNT^ is through the calling arguments. 

3. 3. 4 . 3 Input 

This program reads the ACA index and data files. Control is 
from an interactive terminal. 

3. 3.4.4 Output 

The program writes to the interactive terminal. An optional 
Line Printer report (fig 3. 3. 4. 4) may be generated. 

3. 3. 4.5 Storage 
1683 Bytes 

3. 3. 4. 6 Description 

The ACA file names are requested and the ACA files are opened. 
The optional line printer flag is set, and a segment number is 
requested. If the segment is in the data base, the dates are 
optionally printed. Then a new segment is requested. V7hen 
setment 9999 is entered, the report is optionally printed 
and the program stops. 
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3.3.5 INITIAL PILE GENERATION (SEED) 

3. 3. 5.1 Linkage 

This program calls no subroutines. 

3. 3. 5. 2 Interface 
None 

3. 3. 5. 3 Input 
None 

3. 3. 5 . 4 Output 

Seed creates the ACA files SEEDINX.ACA and SEEDATA.ACA. 

3. 3.5.5 Storage 
923 Bytes 

3. 3. 5.6 Description 

In order for the CDDDS programs to function there must be pre- 
existing ACA files. SEED creates such files with only one data 
record having segment number 1, which is marked for deletion. 
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3.3.6 TABLE OF CONTENTS (TOC) 

3. 3. 6.1 Linkage 

This program calls GETNAM, and also its entry point GETCRD. 

3. 3. 6. 2 Interface 

The interface with GETNAM is through the calling arguments. 

3. 3. 6 . 3 Input 

This program reads the data file and the index file. GETNAM either 
reads FN.DAT, or queries the interactive terminal. 

3. 3. 6 . 4 Output 

The output is a listing of all the files in the ACA data base 
(figure 3. 3. 6. 4). This list is ordered by segment number. 

In debug mode, only the segment numbers are listed. 

3. 3. 6. 5 Storage 
3424 Bytes 

3. 3. 6. 6 Description 

The program first calls GETCRD, an entry of GETNAM, to get the 
index and data file names. Then it opens the index and data 
files. Up to 1024 segment nuiri>ers are read in and sorted. Then 
the data for each segment is read from the data file and written 
to the line printer. 
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3.3.7 TAPE LOAD (TPLOAD) 

3 . 3 . 7 . 1 Linkage 

TPLOAD calls the entry GETCRD of GETNAM and POT. In addition 
several, system routines inclav.od in QIO are called. 

3. 3. 7. 2 Interface 

The interface with GETNAM is through the c lling arguments. 

The interface with PUT is through the calling arguments and 
through the common block/ANCIL/. 

3. 3 . 7 . 3 Input 

One input card is read from the default unit 1 file F0ROO1. 

COL 1-3 Tape drive identifier one of MTO, MTl, XTO. XTl 
COL 5-7 Number of segments to be read in, e.g. 999. 

COL 10-15 Code name to be iiped for loading, e.g SECH**. 

3. 3 . 7 . 4 Output 

Data is transferred from the tape to the .ACA data base. At 
the end, a short line printer report is written (fig 3. 3. 7. 4). 

3 . 3 . 7 . 5 Storage 
5214 Bytes 

3 . 3 . 7 . 6 Description 

TPLOAD reads the data for one acquisition from the tape using 
QIO routines, loads the data into the appropriate array location 
and PUTS the data into the data base. After all requested data 
has been loaded, a line printer report is created. 
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3.3.8 DATA BASE REBUILD (XFER) 

3. 3. 8.1 Linakge 

XFER calls GETNAM's entry point GETCRD. 

3. 3. 8. 2 Interface 

The interface with GETNAM is through calling arguments only. 

3. 3. 8 . 3 Input 

The ACA dat file specified in FN.DAT is read. 

3. 3.8. 4 Output 

All data records in the data file except those marked for 
deletion (byte 254 set .TRUE.) are entered into the new data 
base SEEDINX.ACA, SEEDATA.ACA. 

A small report is generated (fig. 3. 3.8.4) 

3. 3. 8. 5 Storage 
2198 Bytes 

3. 3. 8. 6 Description 

The ACA data file is read directly and sequentially. If the 
segment is marked for deletion, (byte 254 set .TRUE.) the next 
record is read. Otherwise statements copied from PUT are used 
to enter the record into the SEED files. 

If a new smaller data base is required based on factors other 
than duplication, line 28 of the listing is the place to insert 
the logic. 
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3.4 SUBROUTINES FOR LANDSAT DATA MANIPULATION 


3.4.1 SATELLITE CALIBRATION (CALIB (SAT ,NUW) ) 

3. 4. 1.1 Linkage 

Subroutine CALIB does not call any other subroutine. 


3. 4. 1.2 Interface 

The calling arguments are the satellite identifier (SAT) and 
the number of points to be processed, (MUM) . 


4-channel LANDSAT data is passed through the common block /DATA/. 


3. 4 . 1. 3 Input 

N.A. 

3. 4. 1.4 Output 

N.A. 

3. 4. 1.5 Storage 
1025 Bytes 

3. 4. 1.6 Description 

The program recalibrates the contents of the common block /DATA/ 
to the LANDSAT 2 pre-launch calibration. 


Desired recalibration 


Sat identifier 


LANDSAT 2 Post launch -1 
No correction 0 
LANDSAT 1 1 
LANDSAT 2 Pre launch 2 
LANDSAT 3 (Lockheed version) 3 
LANDSAT 3 (I!RIM version) 4 
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3.4.2 DATE CONVERSION (DDATE(YEAR, MONTH , DA Y , YY , DDD , MK ) ) 

3. 4. 2.1 Linkage 

DDATE calls no other program. 


3. 4. 2. 2 Interface 

All information is passed through the calling arguments, e.g. 
YEAR = 1979, MONTH = 'JAN' DAY = 1, YY * 79, DDD * 001. 

3 . 4 . 2 . 3 Input 
None 

3. 4. 2. 4 Output 
None 

3. 4. 2. 5 Storage 
437 Bytes 

3. 4. 2. 6 Description 

For MK = 0, DDATE converts Julian dates to calendar dates, i.e. 
79001 1979, Jan, 1 

For MK = 1 DDATE converts calendar dates to Julian dates, i.e, 
1979, Jan, 1 79001 

The permitted values for MONTH are 'JAN', 'FEB', 'MAR', 'APR', 
'MAY', 'JUN', 'JUL', 'AUG', 'SEP', 'OCT', 'NOV, 'DEC'. 
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3.4.4 DATA ACCESS GET (SEG,ARKAY, CODE) 

3 . 4 . 4 . 1 Linkage 


In case of error, GET calls subrovitine GETNAM through its entry 
GETCRD. 


3. 4. 4. 2 Interface 

GET interfaces with the calling routine through the calling 
arguments, and the common blocks /ANCIL/, /PILNAJ^/, and /DATA/. 

The interface with GETCRD is through the calling arguments. 

3.4.4 . 3 Input 
NA 

3. 4. 4. 4 Output 

When a segment is not in the data base, a line of output is sent 
to the line printer. 


3. 4. 4. 5 


Storage 
2385 Bytes 


3. 4. 4. 6 Description 

GET first opens the ACA files. In case of error opening the 
file, GETNAM is called and it tries again. With a second failure 
it returns. 

After the files are opened, the index file is searched for the 
segment. If the segment is found, the data file for that 
segment is searched for the correct code name. If the sixth 
byte of the code name is the date is also checked. If the 

fifth byte is four records are loaded into the common 

/DATA/. Otherwise the dot data is passed through the argument 
ARRAY. The common block /ANCIL/ is used for ancillary data. 

If data is not found, status is set to 1. If status is already 
1, a message is printed. 
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3.4.5 ACA FILE NAME INPUT (GETNAM (PDN ,PIN) ) 

3. 4. 5.1 Linkage 

GETNAM does not call other subroutines. 

3. 4. 5. 2 Interface 

GETNAM is interactive. The interface with the calling program 
is only through the arguments. 

The entry GETCRD (FAN,FDN,FIN,KQ) is not interactive except when 
an error is encountered. 

3. 4 . 5 . 3 Input 

See figure 3. 4. 5. 4 for an example of interactive file definition. 

3. 4. 5. 4 Output 

See figure 3. 4. 5. 4 for an example of interactive fixe definition. 

*• 

3. 4 . 5 . 5 Storage 
655 Bytes 

3. 4. 5. 6 Description 

GETNAM establishes the ACA data file name FDN and the ACA 
index file PIN. 

Interactively if you enter [a,b] FDN and FIN are defined. 

FDN * [a,b]DOTDATA.ACA. 

FIN = [a,b]DOTINDEX.ACA. 

If you enter FN, file [200,7] FN.DAT is opened and FDN and FIN 
are read. If Z is entered FND and PIN are obtained interactively. 
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For entry GETCRD, the file name in FAN of length KO-1, is 
opened and FDN and FIN are read. If KQ = 1, [200,7] FN.DAT is 

opened and FDN and FIN are read. 

In every case, the arrays PIN and PDN' are zero filled 

before control is returned to the calling program. 
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3.4.6 HAZE CORRECTION (HAZCOR) 

3.4. 6.1 Linkage 
HAZCOR calls KAUTH. 


3. 4. 6. 2 Interface 

HAZCOR interfaces with the calling program through the common 
blocks /ANCIL/, /DATA/, and /KAUTH/. 

The interface with KAUTH is through the common block /KAUTH/. 

3. 4 . 6 . 3 Input 
None 

3. 4 . 6 . 4 Output 
None 

3. 4. 6. 5 Storage 
2245 Bytes 

3. 4. 6. 6 Description 

HAZCOR takes the value of GAM from the common block /ANCIL/, 
computes correction terms and applies them to the contents of 
/DATA/. Then KAUTH transforms the corrected values stored in 
DATA, updating /KAUTH/. 


HAZCOR 

3. 4. 6. 7 Flowchart 
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3.4.7 INTERACTIVE SEGMENT CHOICE (INTRAC (SEG,DATE,CNAM) 

3 . 4 . 7 . I Linkage 


INTRAC calls GETNAM. 

3 . 4 . 7 . 2 Interface 

The interface with the calling program is through the arguments 
and the common block /PILNAM/. 

The interface with GETNAM is through the arguments. 

3 . 4 . 7 . 3 Input 

INTRAC is interactive. See figure 3. 4. 7. 4 for examples. 

3.4 .7.4 Output ; 

INTRAC is interactive. See figure 3. 4.7. 4 for examples. 

There is an optional line printer report. 

3. 4. 7. 5 Storage 
1566 Bytes 

3. 4. 7. 6 Description 

INTRAC makes it possible to interactively choose segment acquisi- 
tions for further processing. 

INTRAC first opens the ACA files. If this fails, it interactively j 

gets the file names through GETNAM and opens the ACA files. It i| 

then asks for a segment number. If it gets 9999 it returns. If 
it gets a segment in the data base, the acquisitions are optionally || 

displayed. If the segment is not in the data base the program j 

asks for a new segment . j 
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Figure 3. 4.7,4 
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3.4.8 KAUTH TRANSFORMATION (KAUTH(N)) 

3. 4. 8.1 Linkage 

KAUTH does not call any other subroutine. 

3. 4. 8. 2 Interface 

KAUTH transforms the first N sets of channel values in the 
common block /DATA/ and puts the result into the common block 


/KAUTH/. 


3. 4. 8. 3 

Input 


None 

3. 4. 8. 4 

Output 


None 

in 

• 

00 

• 

. 

Storage 


1824 Bytes 

3. 4.8.6 

Description 


KAUTH applies the LANDSAT 2 version of the Kauth-Thomas trans- 
formation to the 4 channel data in /DA.TA/ and stores the 4 channel 
result in /KAUTH/. 







Fxgure 3. 4, 8. 8 


Listing for KAUTH 
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3.4.9 MEAN VAIUE COMPUTATION (MEAN (IX, FLAG, N,NN, SEAR, STDEV) ) 

3. 4. 9.1 Linkage 

No other subroutines are called. 

3. 4. 9. 2 Interface 

All values are passed through the calling arguments. 

3 . 4 . 9 . 3 Input 
None 

3. 4. 9.4 Output 
None 

3. 4. 9. 5 Storage 
566 Bytes 

3. 4. 9. 6 Description 

Mean computes the mean and standard deviation of the first N 
values of IX for which the value of FLAG is 1. NN is the number 
of values with flag = 1. XBAR is the mean. STDEV is the 
standard deviation. 
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3.4.10 SPECTRAL PLOT (PLOT(lX,iy,CHAR,N) ) 

3.4.10.1 Linkage 

PLOT does not call any other subroutine. 

3.4.10.2 Interface 

All datt is passed through the calling arguments. 

3.4.10.3 Input 
None 

3.4.10.4 Output 

A spectral plot is written to unit 8 and printed, (figure 3.4.10.4). 

3.4.10.5 Storage 
3215 Bytes 

3.4.10.6 Description 

The first N characters in the BYTE array CHAR are plotted at the 
points (IX, lY) for 20 < 120, -5 < lY < 50. 

Out of range points are truncated and plotted. 
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Listino for PLOT 
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3,4.11 STORE NEW DATA ( PUT ( SEG , ARRAY , NAME ) ) 

3.4.11.1 Linkage 
GETNAM may be called. 


3.4.11.2 Interface 

The interface with the calling program is through the calling 
arguments and the common block /ANCIL/. 

3.4.11.3 Input 
None 

3.4.11.4 Output 
None 

3.4.11.5 Storage 
2022 Bytes 

3.4.11.6 Description 

Put opens the ACA files, then looks for a previous entry for 
segment SEG. If no such entcy has been made, a new entry is made, 
If segment SEG does occur, PUT checks that the current data is 
not duplicated, and if it is not, the current data is added to 
the data base. 


ARRAY is the dot data. NAME is the identifying code name. DATE 
and FILL are passed through the common block ANCIL. 


3-116 








J. 

O; 

X UJ 

^ n* 

C 3 V X 
U If 

c/) > 

It V- •> 

^ ♦ a' 


111 at 9 -^ 

c\ 

•^ /let. 


* lA ,c 
u> 4 i 

' CM < U. 

^ U. 
rt A *:j 
oi - 10 




























L!*CKS/hf> 


m 

in 

lu 










































3.4.12 SCREEN ALL DOTS (SCRALL(N)) 

3.4.12.1 Linkage 

SCRALL does not call any other subroutine 


3.4.12.2 Interface 

SCRALL interfaces with the calling routine through the calling 
argument and the common blocks /ANCIL/, /DATA/, /KAUTH/, /FLAG/. 

3.4.12.3 Input 
None 

3.4.12.4 Output 
None 

3.4.12.5 Storage 
2579 Bytes 

3.4.12.6 Description 

SCRALL is an implementation of the SCREEN algorithm developed for 
LAC IE by ERIM. 

The data are transformed to KAUTH coordinates and stored in 
/KAUTH/. The sun angle correction factor is applied and the 
data are screened. The screening result for each pixel is stored 
in the BYTE array FLAG. A soil line estimate is made, and if 
SLN = 99 it is returned to the calling program. 


SCRALL must be compiled with the CO: 50 switch. 
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3.4.13 DATA SORT (SORT (X , Y,P ,N) ) 

3.4.13.1 Linkage 

Sort does not call any other routine. 

3.4.13.2 Interface 

All data is passed through the calling arguments. 

3.4.13.3 Input 
None 

3.4.13.4 Output 
None 

3.4.13.5 Storage 
159 Bytes 

3.4.13.6 Descriptio 

Sort defines the pointer array P so that X(P(1)) X (P (2) ) <^* • *<X (P (N) ) . 
In case X(P(K)) = X(P(K+1)) then Y(P(K)) < (P(K+1)). 
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3.4.14 TRAJECTORY PLOT (TRPLT (IX , JY ,NT,ND,NB) ) 

3.4.14.1 Linkage 

TRPLT calls no other subroutines. 

3.4.14.2 Interfe^ce 

The interface with the calling routine is through the calling 
arguments and the common block /TRAP/. 

3.4.14.3 Input 
None 

3.4.14.4 Output 

TRPLT creates line printer plots. (fig 6.4.14.4) 

3.4.14.5 Storage 
14,445 Bytes 

3.4.14.6 Description 

TRPLT creates plots of size IX by JY, packed several to the page 
with NB dots per plot, each dot plotted NT times for a total of 
ND (209) dots. The plot format is given by the array PLT. 

Current (Aug 1979) practice is 2 dots per plot, up to 16 times 
per plot, 2 plots per page. The array PLT is read from [200,7] 
GRINT.DAT. 

The points 

(T(I) ,X(I,N)) and 

(T(I) ,Y(I,N) ) are plotted 

with symbol CHAR(I,N). The number of N's per plot is NB. There 
may be several plots per page. 

Changing the array PLT will radically change the output. 
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3.4.15 VEGETATION INDEX (VIN (VEG,ID,N) ) 

3.4.15.1 Linkage 

VIN does not call any other program. 


3.4.15.2 Interface 

The interface with the calling arguments and the common block 
/DATA/. 

3.4.15.3 Input 
None 

3.4.15.4 Output 
None 

3.4.15.5 Storage 
4160 Bytes 

3.4.15.6 Description 

VIN computes the vegetation index specified by ID for the first 
N points of /DATA/ and puts the result in VEG. The identifiers 
are 

1 Greenness 

2 Green Number 

3 Ashburn Vegetation Index 

4 Perpendicular Vegetation Index 

5 Transformed Vegetation Index 6 TVI6 

6 Transformed Vegetation Index 7 TVI7 

7 Kannemasu's simple leaf area index 

8 Kannemasu's complex leaf area index 

9 Berkely ratio 2 *0114/0112 

10 DVI 

11 RVI 

12 Unused 
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3.4.15,7 Flowchart 
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